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Abstract. In this paper we construct the Quantum GauB Jordan Elimination (QGJE) Algorithm 
and estimate the complexity time of computation of Reduced Row Echelon Form (RREF) of an 
N x N matrix using QGJE procedure. The main theorem asserts that QGJE has computation time 
of order 2^/2. 



1. Introduction 
Let us consider a general system of linear equations 

{o-ii^i H h ai n x n = bi 
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we have the system in the matrix form 

(1.3) Ax = b, 

with the augmented matrix A = [A\h]. 

It is well known in linear algebra that by using elementary transformations we can reduce the 
augmented matrix of the system to the so called Reduced Row Echelon Form (RREF) and after 
reenumerate the variables we can rewrite it in the form 



(1.4) 



RREF ( [A | b]) 
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The solutions of the system is some affine m-fiat, m = n — r, r = ikA with a vector basis 
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(1.5) 



and an affine point 



(1.6) 
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In the particular case, where r = n = N the system has a unique solution. Start from here we 
consider the only this nondegenerate case. The procedure of producing the RREF of the augmented 
matrix A = [A\b] is written as some computer program. We refer the readers to [PTVF] for a 
program in C language. It is easy to show that this computation needs at least the time of order 2^. 
The most time complexity is paid to find a pivotal element. 

In this paper we try to simulate this procedure on quantum computers, and show that we could 
have the complexity time of order 2 N I 2 . 

The paper is organized as follows. We state our QGJE Algorithm and the Main Theorem in 
§2. In §3 we analyze the structure of quantum addition and explain its mechanism of work. Its 
seems to the authors that by using the representation theory and the corresponding Fast Fourier 
Transformation is easy to explain the quantum adder and also quantum multiplication. In §5 we 
prepare the main background of quantum computing, and finally in §5 the main theorem is proved. 

2. QGJE Algorithm 

Step 1 Use the Graver's Search algorithm to find out the first non-zero an ^ 0. 

Step 2 If the search is successful, produce the first leading 1 in the first place as an, else change 

to the next column and repeat step 1. 
Step 3 Eliminate all other entries a^i, . . . , a^ i in the column. 

Step 4 Change N to N — 1, control if still N > 0, repeat the procedure from the step 1. 

Step 5 In backward eliminate all ajv-i,Ar, • ■ ■ , cli,n- 

Step 6 Check if iV > 0, change N to N — 1 and repeat the step 5. 

Theorem 2.1. In the Quantum Gaufi-Jordan Elimination Algorithm one needs at most 

N(N - 1)(2JV + 1) 



(2.1) 



+ 



V2-1 



0{2 N ' 2 ) 



operations. 
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Proof of this theorem is the main goal of the paper. 

3. A STRUCTURAL COMMENT OF QUANTUM ADDER CONSTRUCTION 

In this section we give a brief review of quantum computers, as we need in the rest of this paper. 

3.1. g-digit counting system. Let us first remind the ordinary rule of arithmetic operations. Let p 
be a prime and r some positive integer. Write a number in series of q = p r with coefficients in the 
finite field F g , 

m n 

(3.1) a = ^a;g\ b = ^ tyq 3 , a u bj G F q . 

i=0 i=0 

The sum a + b and the product a. b are defined as usually 

(3.2) a + b= (a + a x q + ...) + (b + biq + . . . ) = c + c±q + . . . , 

(3.3) a. b = (a + a±q + . . . )(b + b±q + . . . ) = d + d±q+ : . . . 
where 

(3.4) Co = a + b ( mod q), 

(3.5) Ci = ai + 6i + [(a + b )/q]( mod g), 

(3.6) do = ao.bo(mod q), 

(3.7) di = ai.&i + [(a .6o)/?]( mod g), 



3.2. Qubits. A fundamental notion of quantum computing is the notion of gMfoY. 



Definition 3.1. A qubit is a quantum system the Hilbert space of its quantum states is the q- 
dimensional complex plane C q , and therefore a quantum state of a qubit is a normalized complex 
vector v£C, 



(3.8) 



v = 



v 



Vq-1 



, ||v|| = \v \ 2 + ... + |^_i| 2 = 1. 



Choose the standard basis of q vectors 





1 

An arbitrary vector can be decomposed in this basis as a linear combination 
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(3.9) |0> = 
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(3.10) 



v = u o |0)+ui|l) 
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Lemma 3.2. Every n-digit integer number can be written as a tensor product ofn qubits 

(3.11) a = a <g> • • • <g> a n _! G (C 2 )®",^ eF, = Z/gZ. 
PROOF. Every integer can be written in binary form 

(3.12) a = a + a{l H h a n 2 n , a* G F q = {0, 1, . . . , q - 1}. 

Therefore we have \a ) (g> • • • ® |a n ) G (C 2 ) n . □ 
Let us consider the additive groups G = (F q ) n . A unitary character is a continuous homomor- 
phism x '■ G — > S 1 = U(C) from the group G into the group of unitary automorphism of C, i.e. a 
continuous map \ '■ G — > S 1 = U(C), such that 

(x(x + y) = xWx(y) 
x(o) = 1 
llxWII = i. 

Lemma 3.3. Every (multiplicative) unitary character of the additive groups G = (F q ) n if of the 
form 

2ixi n 

(3.14) Xo(a:) = exp(— ^ a i x i)- 

Q i=i 

PROOF. The lemma is an easy exercise from the group representation theory. □ 
It is convenient to write the number in fractional form, i.e. write 

(3.15) 0.a = ^ + ^i T + .... 

q n q n-l 

The integers modulo 2 n can be written as 

(3.16) k i-> e 2wi0 - k , where O.k = k/q n , 

for k — 0, 1, . . . , q n — 1. Let us consider the Haar measure /j,(k) = -j^. Then we have the natural 
Fast Fourier Transform (FFT). 



Definition 3.4. The transformation 

q n -l 

2 r~ 



q n -l 

(3.17) FFT: |y> - -^e^^k) 



k=0 



is called the Fast Fourier Transform of |y) 

After that we can formally write a number \y) in the form of exp(27ri0.|/ n |/ n _i . . . y ), i.e. we 
have 

(3.18) FFT : \y) ^ exp(2«0.^„_ 1 ...y ) 
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3.3. Quantum addition. By using FFT we can write integers in the form |3.18l i.e. we have a 1-1 
correspondence 



(3.19) 



a FFT ) g27ri0.a n a n _i...ao 



(3.20) 



5 FFT ) e 27ri0.fe„6 n _i...b 



Definition 3.5. The sum of two numbers a and b is defined as follows. 

(3.21) a FFT , e 2ni0.a n ...a add 0-0-0&0, e 2 7 ri0.a n a„_ a ...a 1 (a +&o) add 0-0-Ofri, 



(3.22) 



add e 27riO.(a n +b n )...(ao+f> ) 



FFT inverse 



a + b 



Remark 3.6. The reduction (memorize [(a* + bi)/q] to a i+ i + b i+ i is the phase transition with the 
matrix 



(3.23) 



1 

e 



2tti 







3.4. Quantum gates. Let us recall some fundamental gates in quantum computing. 



3.4.1. Hadamard gate. 



Definition 3.7. The 1-qubit Hadamard gate — H — is defined by the matrix 
(3.24) H 



s/2 



1 1 

1 -1 



The Hadamard gates acts on 1-qubits as follows 
(3.25) 



\x) 



H 



3.4.2. Phase gate. 



Definition 3.8. The phase 2-qubit gate — <J> — is defined by the matrix 



(3.26) 

it produces an action 
(3.27) 



V2 
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e J * 
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— e % ®\x). 
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3.4.3. CNOT(XOR) gate. 

Definition 3.9. The 2-qubit CNOT (XOR) gate is defined by the matrix 

10 



(3.28) 

The action of this gate is given by 
(3.29) 

where x © y := x + y mod (2). 
3.4.4. Unitary gates. 



C 



\x) 

\y) 



10 
1 
10 



\x) 

\x@y), 



Definition 3.10. The unitary gate - U - is defined by a unitary matrix U, 
(3.30) 



U 



1 

o u 



3.4.5. Problem of quantum computing. One of the major result concerning this problem is the 
following classification result. 

Theorem 3.11 (Universal gates). One-qubit gates, CNOT gate, phase gate and one of unitary gate 
generate all the other gates. 

4. Quantum Grover's Search Algorithm 

4.0. 6. Black-box ( Query) problem. The problem is to compute some Boolean function of type 
(4.1) /:{0,l} n ^{0,l}. 

4.1. Deutsch's original problem. The problem is to check whether the above defined Boolean 
function is balance: the number of values 0's equals to the numbers of values l's, or constant. 



4.2. Search problem. Given a Boolean function / : {0, 1}™ — > {0, 1}, defined by fk(x) 
Find k. 



Jxk- 



4.3. Deutsch's scheme — |"h] — \t\ — [h]— . The Deutsch's scheme produces the effects (modulo 
a constant, we often omit it in order to not the formula longer, but it is easily recovered in the case 
of necessary): 

(4.2) |*>(|0) - |1» ^ U ' (-lY^l (|0> - |1». 

measurement 

In this scheme, the second qubit is of no interest, while the first qubit has state 



(4.3) 



(_l)/(o)|o) + (_i)/(D|i) = 



±(|0) + |1)), if / = const 
±(|0) — jl)), if / = balanced 
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4.4. Quantum Grover's Search Algorithm. Scheme: — f k 
This algorithm has special effect: Choose the state 

1 



H 



/o 



H 



(4.4) 

and produce the Grover's iterate 
(4.5) 



i=i 



I^>(|0>-|1» 



H 



/o 



H 



| measurement) ( 1 0) — |1)). 



This iterate does nothing to any basis element \i) except for \k) is changing into —\k), i.e. the 
reflection in the hyperplane perpendicular to the plane 



(4.6) 



N-l 

H|0> = |5> = ^El < >. 



i.e. rotation in the plane generate by \k) and \S) an the Grover's iterate is a rotation of twice the 
angle from {k 1 ) to S 1 ). 

5. Proof of the Main Theorem 

an 



Step 1. Let us consider the first column 



. Each an is either zero or non-zero. We 



. ClNl . 

do find the first non-zero value. In classical computation we needs to check all possible values of 
column entries and we needs 2 N operations. In quantum Grover's Search Algorithm we need to 
proceed 2 Ar / 2 operations. More precisely, choose a state 



JV-1 



i=0 



(5.1) |S) = 2-"/ 2 £ l<> 

and produce the Grover's iterate 
(5.2) 



l*>(|o>-|i» 



fk 



H 



/o 



H 



measurement) ( 1 0) — |1)). 



This iteration does nothing to any basis element except for \k) is changing into —\k), i.e. reflection 
in the hyperplane perpendicular to plane 

(5.3) H |0) = 2- Ar / 2 ^K) = |5), 

i=0 

i.e. a rotation in the plane generated by \k) and \S), and the desired state \k) exactly is $, the 
iterate should be repeated m times with 



(5.4) 

Since 

(5.5) 



/~ , x , 71 7T 1 

2m + 1 $ — or m w — . 

v ; 2 4$ 2 



—^77- = sin $ w $, we have m — 2 N ^ 2 . 

2N/2 4 
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This estimate is optimal in the sense that any another quantum algorithm for searching an unstruc- 
tured database must stake time at least of order 0(2 N I 2 ). 

Suppose we have found a non-zero element and denote it (the first nonzero element) by an 

Step 2. Produce the first modified row 

(5.6) ri = [1, ]. 

an an 

For this we need at least N — 1 multiplications with Register the resulted row as the first row. 
Step 3. 

• Multiply the new first row with — aji. For this we need N — 1 operations of multiplication 
and 

• Add the result to the j th row. For this we need another N — 1 operations. 

Totally in this step we need 2(N — l) 2 arithmetic operations. 

Step 4. Now change N to N — 1 (need 1 subtraction) and control for the value N > 0. We 
need one CNOT gate. 

Loop back to the first step produce summation from 1 to N. 

Steps 5 and 6. Eliminating backward the numbers cln-i,n, • • • , a>i,N, we need iV— 1 additions 
and one CNOT gate to control and change N into N — 1 . 
All Steps. Totally in all steps we need at most 

N 

(5.7) ^(2 n/2 + 2(n-l) 2 + 2(n-l) + l + l) 

n=l 

operations. 

Reduction. The following lemma is trivial 
Lemma 5.1 (Combinatorial Lemma). 



(5.8) J2 



2 N(N + l)(2N + l) 



n 



6 

Finish the Proof of the Main Theorem Denote the integral part of a real number a by 
[a] . Finally we have 

N 

Y^(2 n/2 + 2(n - l) 2 + 2{n - 1) + 1 + 1) w 

n=l 

N N N 

2 E n+2 E( n - 1 ) 2 + E v/ 2 n ] = 

n=l n=l n=l 

2 N(N±21 + 2 N-1)N(2N- 1} ^ V2"-l = 
2 6 1 V2-1 

(jV - 1)7V 3 (2jV+1) + [V2^1 » 0(2^ 

Corollary 5.2. For N >> 1, the number of operations needed to produce the RREF of a square 
matrix using QGJE Algorithm is of order 0(2 N / 2 ). 
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In the next paper HDGI we apply this computation to simulating of accounting principles on 
quantum computers. 
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